Proposal of a robust measurement scheme for the non-adiabatic spin torque using the displacement 

of magnetic vortices 
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The strength of the non-adiabatic spin torque is currently under strong debate, as its value differs by orders 
of magnitude as well in theoretical predictions as in measurements. Here, a measurement scheme is presented 
that allows to determine the strength of the non-adiabatic spin torque accurately and directly. Analytical and 
numerical calculations show that the scheme allows to separate the displacement due to the Oersted field and is 
robust against uncertainties of the exact current direction. 
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A spin-polarized current flowing through a ferromagnetic 
sample interacts with the magnetization and exerts a torque on 
the local magnetic moments. For conduction electron spins 
that follow the local magnetization adiabatically it has been 
shown that the interaction via spin transfer can be described 
by adding a current-dependent term to the Landau-Lifshitz- 
Gilbert equation. HI This equation has been extended by an 
additional term that takes the non-adiabatic influence of the 
itinerant spins into account. Q Theoretically, several mecha- 
nisms have been proposed as the origin of the non-adiabatic 
spin torque, leading to different orders of magnitude for its 
strength. Q IH 0, H 01 Thus a precise measurement of the 
non-adiabatic spin torque is necessary to give insight into its 
microscopic origin. A determination of its strength is fur- 
ther important for a reliable prediction of the current-driven 
domain-wall velocity. 101 Currently measured values of the 
non-adiabatic spin torq ue for permalloy differ by one order of 
magnitude. ll7[ Isl l9[ llOll thus the strength of the non-adiabatic 
spin torque is under strong debate. In these experiments the 
observed motion of a domain wall was compared with mi- 
cromagnetic simulations to determine the non-adiabatic spin 
torque. This analysis is highly susceptible to surface rough- 
ness and Oersted fields. 

Due to its high symmetry and spacial confinement a vor- 
tex in a micro- or nanostructured magnetic thin-film element 
is a promi sing sy stem for the investigation of the spin-torque 
effect, i 111 I I2II 1311 Vortices are formed when the in-plane mag- 
netization curls around a center region. In this few-nanometer- 
large center region, called the vortex core, the magnetization 
turns out of plane to minimize the exchange energy. There 
are four different ground states of a vortex. These states are 
labeled by the direction of the out-of-plane magnetization, 
called polarization p, and the sense of rotation of the in-plane 
magnetization, called chirality c. Polarizations of p = 1 and 
p = — 1 denote a core that points parallel or antiparallel to the 
z axis, respectively. A chirality of c = 1 denotes a counter- 
clockwise curling of the in-plane magnetization while c — —I 
denotes a clockwise curling. 

It is known that vortices are displaced from their equilib- 
rium position when excited by spin-polarized electric current 



pulses. iH [llHIilEIlTlIillillllThe spatial con- 
finement of the vortex core within the film element yields 
an especially accessible system for measurements with scan- 
ning probe techniques, such as soft x-ray microscopy, x-ray 
photoemission electron microscopy, or scanning electron mi- 
croscopy with polarization analysis. An analytical solution of 
the extended Landau-Lifshitz-Gilbert equation shows that for 
a current-driven vortex the forces due to the adiabatic and the 



non-adiabatic spin torque are perpendicular to each other 11 15 1 
In this work we present a scheme which allows to measure 
the contributions due to the adiabatic spin torque, the non- 
adiabatic spin torque, and the Oersted field separately. It bases 
upon analytical calculations lIlSll and overcomes the two main 
difficulties that occur in an experiment. The first problem 
arises from an additional vortex displacement due to the Oer- 
sted field accompanying the current flow. 1211 This displace- 
ment is comparable in size to the displacement due to the non- 
adiabatic spin torque and both displacements point in the same 
direction. (1151 Thus, the unknown contribution of the Oersted 
field has to be separated from the measured signal. The sec- 
ond problem is the exact determination of the displacement 
angle. Since the displacement due to the adiabatic spin torque 
is about one order of magnitude larger than the displacement 
due to the non-adiabatic spin torque, a small uncertainty in the 
direction of the current through the sample would cause large 
errors in the determination of the non-adiabatic spin torque 
parameter. To test the applicability of our analytical findings 
they are applied to vortex displacements obtained from three- 
dimensional micromagnetic simulations. 

For the analytical calculations we start from a modified ver- 
sion of the Thiele equation 12 ll |22 1 



F + Goe, X (v, + bjj) + Drav^ + Do^bJ = (1) 

that takes deformation of the vortex into account, jl^ Here 
Vc is the velocity of the vortex core, a is the Gilbert damp- 
ing, ^ the strength of the non-adiabaticity, F the force on the 
vortex. Go the z component of the gyrovector, and Dq the 
diagonal element of the dissipation tensor. The coupling con- 
stant bj = Pfis/ieMs) between the current and the magne- 
tization depends on the saturation magnetization and the 
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spin polarization P. The assumption of a magnetization pat- 
tern which rigidly gyrates holds true only for the small vortex 
core. Due to the spatial confinement the remaining part of the 
vortex has to deform while the core is moving. Dp < Do is a 
phenomenological parameter that takes into account a reduced 
dissipation due to this deformation. 1 23] 

We will investigate a square thin-film element with a ho- 
mogeneous current flowing in x direction. The accompany- 
ing Oersted field is accounted for by a homogeneous field 
H in y direction. For small displacements of the vortex core 
from its equilibrium position the demagnetization energy can 
be expanded up to second order in the core displacement 



R = {X, Y). The force on the vortex is then given by BlSll 



pLoAIsHldc + muj'^X 



(2) 



with the lateral extension I, and thickness d of the system. The 
factor mcjj parameterizes the confining potential. 11 511 

This equation can be solved for harmonic excitations of the 



form H{t) = Hoe By and = joe e,-,. After all tran- 
sient states are damped out, the position of the vortex core is 
given by 
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with H = ^Hol / {2tt), the gyromagnetic ratio 7, j = bj jo, 
the free oscillation frequency of the vortex uj, and its damping 

From Eq. (|3]l it is obvious that an Oersted field has the same 
influence on the vortex as the non-adiabatic spin torque. Thus 
the presence of an Oersted field can disturb the measurement 
of the non-adiabatic spin torque. In experiments the coordi- 
nate system is given by the sample axis. A small uncertainty 
of the direction of the current flow, e.g. due to a rotation or im- 
perfections of the sample, yields a mixing of the displacement 
components, resulting from the adiabatic spin torque and the 
smaller non-adiabatic spin torque, relative to the sample axis. 
This mixing causes a large error in the measurement of the 
displacement originating from the non-adiabatic spin torque. 

An excitation with a direct current causes a displacement of 
the vortex core to a new steady-state position. l24ll A benefit 
is that a direct current allows for a measurement with a non- 
time-resolving technique. The displacement of a vortex driven 
by a direct current is given by 
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The sign of the displacement induced by the Oersted field 
depends on the chirality of the vortex, while the displace- 
ment due to the adiabatic spin torque is determined by the 
polarization. ||20| The non-adiabatic spin torque causes a dis- 
placement that is independent of the vortex properties p and c. 
Vortices with different p and c values can be achieved by re- 
magnetizing the sample. Comparing the displacement of three 
vortices with different polarizations and chiralities it is there- 
fore possible to separate the contributions of all three forces 



to the displacement of the vortex. From Eq. (|4]i we find 
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These equations are schematically illustrated in Fig.[T] From 
Eqs. dSal i and ( Bbl i it is possible to determine the non- 
adiabaticity parameter as 
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Since this equation is independent of the strength of the Oer- 
sted field, the angle of the sample, the frequency ui, the damp- 
ing r, and the parameter Dp, it yields the sought measurement 
scheme. With this scheme a direct determination of ^ is acces- 
sible. Only one micromagnetic simulation for the determina- 
tion of |Do/Go| is necessary since \Do/Go \ is independent of 
£, and j. 

Micromagnetic simulations of the experimental setup al- 
low to determine the positions of the vortex core with a pre- 
cise knowledge of the micromagnetic parameters of the sys- 
tem. The simulations therefore allow to test the analyti- 
cal results in Eqs. (|4|i and For the simulations the ma- 
terial parameters of permalloy, i.e., a saturation magnetiza- 
tion of A/s = 8 • 10^ A/m and an exchange constant of 
A = 1.3 • 10^^^ J/m, are used. Since we are interested only 
in the steady final position of the vortex we used a Gilbert 
damping of a — 0.5 to ensure a fast damping of the transient 
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FIG. 1: Scheme of the determination of the three different contribu- 
tions to the vortex displacement. By measuring the distance between 
the positions of two different vortices it is possible to separate the 
displacements (a) due to the non-adiabatic spin torque, (b) the adi- 
abatic spin torque, and (c) the Oersted field. Points and crosses de- 
note cores with positive and negative polarization, respectively. The 
in-plane magnetization is denoted by the solid arrows. The dashed 
arrows denote the current direction. For the sake of illustration the 
displacements are exaggerated. 



States to reduce computation time. As a sample system we 
considered a square thin-film element of length I = 500 nm 
and thickness d = 10 nm with a cell size of 2 nm in the lateral 
directions and 10 nm perpendicular to the film. This system 
allows for a reasonable computation time. For the simulations 
we used our extended version of the Object Oriented Micro- 
magnetic Framework. li25i, \26A 

Figure 12] shows the displacement of the vortex core in sim- 
ulations without Oersted field. As predicted by Eq. (|4|i the 
displacement in the direction of the current flow is propor- 
tional to ^ and the displacement perpendicular to the current 
flow is independent of ^. From these simulations the value 
\Do/Go\ ~ 2.26 can be determined. 

In experimental samples we are faced with an Oersted field 
and an uncertainty of the direction of the current flow. To 
mimic the Oersted field in the simulations we applied an in- 
plane field perpendicular to the current. The strength of the 
field is given by H/ [jP) = 1 nm. The uncertainty of the 
direction of the current flow was taken into account by rotating 
the sample by 5 degrees. Figure [3j a) shows the positions of 
the vortex core for both simulations. It becomes visible that 
the Oersted field and the rotation of the sample strongly shift 
the core positions, complicating the determination of ^. 

To test the analytical model we compared the non-adiabatic 
spin-torque parameter that was inserted into the Simula- 
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FIG. 2: Numerically calculated displacement of the vortex core due 
to a direct spin-polarized current of density jP in the absence of an 
Oersted field, (a) The displacement parallel to the current is pro- 
portional to ^. (b) The displacement perpendicular to the current is 
independent of ^. The lines are fits with the linear model in Eq. l|4j- 
For large current densities small non-linear effects can be seen. 



tions with the value ^out that was calculated from Eq. ^ using 
the core positions. Here it is worth noting that the value of the 
Oersted field and the angle of the sample are not needed for 
the calculation of ^out- The results are shown in Fig.[3jb). It 
can be seen that aU the perturbations that are inserted in the 
simulations can be effectively excluded by the analytical cal- 
culations. 

In experimental samples we are also faced with the 
anisotropic magnetoresistance (AMR) effect that leads to in- 
homogeneous current paths, i.e., a higher current density in 
the vortex core. Simulations including these inhomogeneous 
current paths yield a small shift to lower values of ^out- This 
shift is up to 2 % for an AMR ratio of 10 %. 

In the remaining part we will discuss the experimental ac- 
curacy in the determination of ^ that can be achieved with the 
presented scheme. In experiments direct currents of densities 
up to 1.5 • 10^^ A/m^ have been realized in permalloy on a 
diamond substrate. 1271] Assuming a spin polarization of 0.5 
we get a spin-polarized current density of 0.75 • 10^^ A/m^, 
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FIG. 3: (a) Position of the vortex core displaced by a spin-polarized 
direct current of density jP = 3 ■ 10^^ A/m with a non-adiabatic 
spin-torque parameter of ^ — 0.1. The overlapping open symbols 
denote the positions for a current in exact x direction without Oer- 
sted field. The closed symbols denote the positions with an applied 
Oersted field and a rotation of the sample by 5 degrees around its 
midpoint (plus). For the latter case the direction of the current is de- 
noted by the arrow, (b) Results for ^out derived from the positions 
shown in (a) using Eq. l[6} for different current densities, ^j,, is the 
value of the non-adiabaticity parameter that was used for the simula- 
tions. 



i.e., the maximum shown in Fig.|2] This yields values of up to 
j = 55 m/s. 

The displacements of the vortex in the numerically inves- 
tigated samples are small compared to the experimental res- 
olutions available. A larger displacement of the vortex can 
be achieved by increasing the lateral size of the structure. 
For example simulations of a square thin-film element of 
length I = 5000 nm and thickness d = 10 nm yielded val- 
ues of \Do/Go\ = 3.82 and uj/{lu^ + T^) = 1 ■ IQ-* s. 
With these values Eq. dSbl i yields 2i?ad =1100 nm. We as- 
sume that the core position can be measured with a resolu- 
tion of (5(2i?non-ad) = 20 nm. Equation (|6]l then yields that 
6^ ~ 0.005 can be realized. This resolution ranges from 5 % 
to 50 % depending on the value of ^.iH i, 11, The resolu- 
tion can be further increased by using thin-film elements with 



still larger lateral sizes. 

In conclusion we present a robust and direct measurement 
scheme for the non-adiabatic spin torque using the displace- 
ment of magnetic vortices. This scheme allows to distinguish 
between the displacements of the vortex core due to the non- 
adiabatic spin torque, the adiabatic spin torque, and the Oer- 
sted field, independently of the direction of the current flow. 
We also showed that an inhomogeneous current due to the 
AMR effect can be neglected. The scheme thus allows a pre- 
cise measurement of the non-adiabatic spin-torque parameter 
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SUPPORTING MATERIAL 

For an analytical investigation the motion of the vortex 
is commonly described employing the Thiele equation. Hsl 
il """" """" 



17, 18 



M escribed employmg the Ihiele equation. yjl 
211 uSL I29I1 This equation is exact for the 
steady state motion of a non-deformable magnetization pat- 
tern. However, this assumption holds true only for the small 
vortex core. Due to the spacial restriction the vortex has to 
deform while the core is moving. This yields a small modifi- 
cation of the Thiele equation that is especially important for 
current-driven vortex motion in view of the non-adiabatic spin 
torque. 

Here we present a modified Thiele equation that takes a de- 
formation of the outer part of the vortex into account. 

With the magnetization M and the magnetic field H a gen- 
eral version of the Thiele equation reads ll2lll 
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(7) 



with the saturation magnetization A/5, the gyromagnetic ra- 
tio 7, the current density j, the Gilbert damping a, the non- 
adiabaticity parameter ^, and the coupling constant hj be- 
tween current and magnetization. 6 and </i are the out-of -plane 
and in-plane angle of the magnetization, respectively. The ve- 
locity V = v{r) of the magnetization pattern may depend on 
the position. Assuming that the magnetization pattern does 
not deform the velocity is independent of the position. Then 
Eq. (|7]) can be written in its well known form|l22 1 



F + G X (uc + bjj) + D{avc + ihjj) 
with the velocity Vc of the vortex core. Here 
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denotes the force on the magnetization pattern. 
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(10) 



is the gyrovector and 



D = -- 



7 



dv {veve + siv?{e)V(jN(j)) (ii) 



is the diagonal dissipation tensor with D^x = Dyy = Do and 
Dzz = 0. The term Davc describes the dissipation of energy 
due to the changing magnetization. 

The integrand in the gyrovector is non-zero only in the 
small vortex core where the out-of -plane angle 9 varies while 
the integrand in the dissipation tensor is also non-zero outside 
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FIG. 4: Scheme of the magnetization (dashed arrows) in a square 
magnetic thin-film element with a vortex. The solid arrows denotes 
the velocity v of the vortex core and of different points within the 
domain wall. 
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FIG. 5: Values of the strength Dr of the dissipation (open symbols) 
and the strength Do of the non-adiabatic spin torque (closed sym- 
bols). The data for films of 10 nm, 20 nm, and 30 nm is denoted by 
squares, circles, and triangles, respectively. 



the core. Close to the boundaries of the sample the magneti- 
zation pattern moves slower compared to the center as it can 
be seen in Fig.|4] Thus the velocity in the third term of Eq. ^ 
depends on space. Aiming a similar form as in Eq. ^ we 
replace the spatially dependent velocity ii in the third term of 
Eq. (|7]i by an effective value Ve which is independent of the 
position. This effective velocity occurs only in the third term 
as the second term is located at the vortex core. For a homo- 
geneous current flow bjj is constant over the sample. Thus we 
do not replace the current by an effective value. The equation 
then reads 



(12) 



The effective velocity Ve depends on the core position R = 
{X, Y) and the core velocity Vc- For small deflections of the 
vortex core, i.e., small deformations of the vortex, Vg can be 
expanded in R and Vc- For Vc = the magnetization is static 
and We = 0. Thus the first non- vanishing term in the expansion 
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FIG. 6: Calculated position of a vortex core excited with a spin- 
polarized direct current of density jP — 6 ■ 10^^ A/m^ in a 1000 nm 
X 1000 nm square thin-film element. Shown is the x position (a) 
and the y position (b) versus time. A film thickness of 30 nm and 
5 = Q = 0.1 was used. The solid (red) line is the vortex core 
position extracted from simulations. The dashed (blue) line is a fit 
with the theory based on the original Thiele equation. 



is proportional to Vc. Here and hereafter we write 



(13) 



Since the effective velocity is always smaller than the ve- 
locity Vc of the vortex core, the new constant Dr is smaller 
than Dq. Inserting Eq. (fT3] l in Eq. (fT2] i yields a modified 
Thiele equation 



F + Gx{vc + bjj) + Dravc + Da^bjj = 0. 



(14) 



Employing the same conversions as used for the original 
Thiele equation 11511 we find an expression for the velocity of 
the vortex core 

{Gl + Dla^)vc = G X F - DtuF - (G^ + DrDaaC)bJ 
+ bj^DoG X j — bjaDrG x j. 

(15) 
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FIG. 7: Calculated position of a vortex core. The solid (red) line is 
the same as shown in Fig. [6] The dashed (blue) line is a fit with the 
theory based on the modified Thiele equation. 



We investigate a square thin-film element with a current in x 
and a field in y direction. With the stray-field energy for small 
deflectionslflal 



and the total Zeeman energy 11511 

Ez = fioMsHldcX 

we get a force of 



(16) 



(17) 



(18) 



Here I and d are the lateral extension and the thickness of the 
square, respectively. As for the original Thiele equation in the 
absence of current and field the excited vortex performs an 
exponentially damped spiral rotation around its equilibrium 
position. The free frequency 



pGomuj^ 
Gl + D^ 



(19) 
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and the damping constant by the frequency and the damping constant. The velocity of 

the vortex then reads 



Gl + Dla^ 



are sHghtly changed compared to their values derived from the 
homogeneous Thiele equation. In the following we express 

Dr=^-^ (21) 
I 
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(22) 



This equation can be solved for harmonic excitations of the form H{t) = H^e Cy and i{t) = jqc 
vortex motion is then given bv lllSll 



eV. The solution for the 
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with H ~ '-fHQl/{2TT) and j = bjjo. The first two terms with 
prefactors A and B are exponentially damped and depend on 
the starting configuration. 

The values of Dp and Dq can be determined by micro- 
magnetic simulations. For these simulations we used our ex- 
tended version of the Object Oriented Micromagnetic Frame- 
work (OOMMF) that includes the adiabatic and non-adiabatic 
spin torque. Hsl The position of the vortex core was 
defined as the point with the maximum out-of -plane magneti- 
zation. To determine this maximum, the simulation cell with 
maximum out-of-plane magnetization and its next neighbors 
are interpolated with a polynomial of second order For the 
simulations the material parameters of permalloy, i.e., a satu- 
ration magnetization of Alg = 8 • 10^ A/m and an exchange 
constant of ^ = 1.3 • 10"^^ J/m where used. 

For the determination of Dr the vortex was excited by a 
magnetic field pulse. The subsequent oscillation was then fit- 



ted with the first two terms in Eq. ( |23] ). Dr can then be deter- 
mined from Eq. (|2TI) . Finally the value of Dq was determined 
by fitting an excitation with a direct current. The results are 
shown in Fig.|5]for different edge lengths I and different thick- 
nesses of the sample. It can be clearly seen that Dr is smaller 
than Dq. Using the original Thiele equation for the descrip- 
tion of the vortex motion Dp and Dq are assumed to be equal. 

Figures|6]and|2]show an example for the fit of a numerically 
calculated vortex-core trajectory using both theories. The the- 
ory based on the modified Thiele equation shows better accor- 
dance than the theory based on the original Thiele equation. It 
can be seen that the Thiele equation has to be modified for a 
sufficient description of the dynamics of current-driven mag- 
netic vortices in the presence of a non-adiabatic spin torque. 
This modification takes the deformation of the outer part of 
the vortex into account. 



